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70.  abstract  (Continue  on  rcvmrmc  cldc  11  nmcmccsry  mnd  Idontlly  by  block  number) 

•“^The  effects  of  surface  roughness  on  the  transition  of  a -laminar  boundary  layer 
are  basically  two-fold:  i)  it  can  distort  the  mean  velocity  and  temperature 
profiles  thus  altering  the  stability  characteristics  of  the  layer  and  ii)  it  can 
induce  additional  disturbances  to  the  boundary  layer  which  may  lead  to  an  earlier 
transition.  Only  the  first  of  these  effects  is  considered  in  this  report.  Three 
different  types  of  surface  roughness  are  investigated  here;  distributed  roughness; 
surface  waviness,  and  a single  two-dimensional  bump.  An  analytical  model  which 
describes  the  effect  of  distributed  roughness  on  the  mean  flow  and  its  stability  - 
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is  described.  This  model  is  further  validated  by  comparison  with  Achenbach's 
experimental  results  and  is  then  applied  to  the  flow  about  an  underwater  body  in 
order  to  illustrate  the  effects  of  surface  roughness  in  a more  realistic  situation 
It  is  found  that  surface  heating  in  the  presence  of  surface  roughness  is  far  less 
effective  at  prolonging  laminar  flow  than  for  the  smooth  wall  situation.  The 
incorporation  of  this  model  into  the  existing  TAPS  code  is  also  described.  For 
surface  waviness  (resulting  from  manufacturing,  perhaps),  Lessen's  model  is 
extended  here  to  provide  a qualitative  assessment  of  its  effect.  The  preliminary 
results  suggest  a weak  effect  for  expected  surface  waviness  characteristics. 

The  third  type  of  surface  irregularity  - a single  two-dimensional  bump  - is 
addressed  using  Smith's  triple-deck  formulation.  Preliminary  analysis  is  pre- 
sented to  provide  the  foundation  for  a future  quantitative  assessment :of  the 
effects  of  surface  roughness  of  this  type,  if  needed.  In  summary>^sur?ace 
irregularities  can  play  a dominant  role  in  determining  the  practical  limit  for 
the  application  of  boundary- layer  control  techniques.  The  present  distributed 
roughness  model  as  incorporated  into  the  TAPS  code  can  provide  a reasonable 
prediction  of  its  effect  in  the  presence  of  boundary-layer  control  techniques 
(shaping,  heating,  suction)  on  the  transition  location  for  an  underwater 
vehicle./*  As  a final  caution,  extreme  care  should  be  taken  in  the  interpre- 
tation ot  any  experimental  results,  if  no  quantitative  surface  finish 
information  is  known. 
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Complex  wave  number  = + 1 

Falkner-Skan  pressure  gradient  parameter 
Parameter  in  roughness  model 
Boundary-layer  thickness 
Displacement  thickness 
Thermal  diffusivity  due  to  roughness 
Momentum  diffusivity  due  to  roughness 

"-■^e 

Wavelength  of  surface  wave 
Normal  coordinate  in  Levy-Lees  transformation 
Radius  of  curvature 
Momentum  thickness 

Angular  location  on  cylinder,  degrees 
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Total /effective  molecular  viscosity 
Laminar  kinematic  viscosity 

Streamwise  coordinate  in  Levy-Lees  transformation 
Density 

Real  frequency 


Subscripts 

e Quantity  evaluated  at  the  boundary-layer  edge 

k Quantity  evaluated  at  the  roughness  height 

w Quantity  evaluated  at  the  wall 


1.0  INTRODUCTION 


Recent  technical  advances  have  successfully  demonstrated  the  feasibility  of  main- 
taining a laminar  boundary  layer  over  the  surface  of  an  underwater  body  operating 
at  a moderately  high  Reynolds  number  (or  speed).  Various  means  of  boundary- 
layer  control  techniques  which  allow  this  capability  to  be  extended  to  substan- 
tially higher  Reynolds  numbers  have  been  developed  and  tested,  or  are  being 
tested,  both  in  small-scale  laboratory  and  in  field  experiments.  The  current 
techniques  for  delaying  the  onset  of  boundary-layer  transition  which,  moveover, 
reduce  the  total  viscous  drag  at  high  Reynolds  number  Include  body  shaping, 
surface  heating  and  suction  through  the  surface.  Depending  on  the  specific 
application,  one  or  a combination  of  these  techniques  can  be  used.  The  design 
of  such  low-drag,  advanced,  underwater  vehicles  relies  critically  on  the  ability 
to  predict  quantitatively  the  various  phenomena  which  might  affect  the  transition 
location.  For  years,  the  most  successful  schemes  for  predicting  the  location  of 
boundary- layer  transition  have  either  relied  on  empirical  data  correlations 
(which  work  well  for  a restricted  range  of  parameters)  or  on  the  semi -empirical 
«en"  criterion  (which  is  based  on  a coupling  of  linear  stability  calculations 
and  empirical  results).  Both  methods  suffer  from  empiricism  which  restricts 
their  validity  to  a limited  parameter  domain.  The  "e”"  method,  in  particular, 
has  been  considered  as  adequate  for  estimating  the  transition  location  on  a 
smooth  surface  over  the  range  of  parameters  for  which  data  exist.  The 
validation  of  the  "e”"  method  for  higher  Reynolds  number  flows  using  various 
boundary-layer  control  techniques  has  not  been  conclusive.  Various  investiga- 
tions, both  analytical  and  experimental,  are  underway  to  provide  a more  reliable 
means  of  predicting  the  location  of  boundary-layer  transition.  At  the  present 
time,  the  engineering  design  tool  most  conmonly  used  is  the  TAPS  program 
(Transition  Analysis  Program  System)  developed  by  McDonnell  Douglas  (Gentry, 
1976).  This  computer  program  does  not  include  the  effects  of  several 
potentially  important  factors  which  are  present  in  a realistic  environment, 
however.  For  example,  the  effects  of  freestream  disturbances  (either  acoustic 
waves  or  turbulence)  and  surface  roughness,  (either  distributed  or  isolated) 
have  not  been  incorporated  into  TAPS  despite  the  fact  that  they,  conceivably, 
could  be  the  dominant  factors  in  many  practical  situations.  In  this  report,  an 


analytical  investigation  of  the  effects  of  surface  roughness  on  boundary-layer 
transition  is  presented. 
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The  effects  of  surface  roughness  on  the  transition  process  are  basically  two-fold: 
i)  Surface  roughness  can  distort  the  mean  velocity  and  temperature  pro- 
files in  a laminar  boundary  layer  from  their  smooth-wall  shapes,  thus 
altering  the  stability  characteristics  of  the  boundary  layer,  and 
ii)  Surface  roughness  can  introduce  additional  disturbances  within  the 
boundary  layer  which  may  lead  to  earlier  transition. 

The  exact  manner  in  which  surface  roughness  affects  the  transition  location  depends 
on  the  type  of  surface  roughness  present.  In  this  report,  three  different  types  of 
surface  roughness  are  investigated;  distributed  roughness,  surface  waviness,  and  a 
single  isolated  "bump." 


In  Section  2,  the  effects  of  distributed  roughness  are  discussed  in  some  detail 
starting  with  a brief  reconstruction  of  the  physical  model.  The  perturbation  to 
the  mean  flow  profile  is  modeled  by  an  enhanced  momentum  and  heat  transfer  near 
the  surface,  which  may  be  considered  as  resulting  from  the  unsteadiness  in  the 
flow  over  the  surface  irregularities.  This  augmented  momentum  and  heat  transfer, 
when  interpreted  in  a statistical  sense,  gives  rise  to  the  concept  of  a "roughness 
sublayer"  situated  between  the  surface  and  the  outer  laminar  portion  of  the  bound- 
ary layer.  This  roughness  sublayer  alters  the  mean  velocity  and  temperature  pro- 
files within  the  laminar  portion  of  the  boundary  layer.  The  distributed  roughness 
model  is  validated  by  comparing  with  Achenbach's  experimental  results  in  Sub- 
section 2.2.  The  model  is  then  applied  to  a realistic  body  shape  to  illustrate 
its  predictive  capability  in  a more  quantitative  manner.  The  incorporation  of 
this  model  into  the  existing  TAPS  code  in  order  to  provide  the  technical  community 
with  the  additional  capability  of  including  vehicle  surface  roughness  conditions 
into  performance  predictions  is  described  in  Section  3. 


In  addition  to  distributed  roughness,  other  types  of  surface  roughness  which  are 
considered  here  are  two-dimensional  waviness  and  a single,  isolated  bump. 
Although  waviness  and  single  roughness  elements  will  affect  transition  through 
the  same  two  basic  mechanisms  mentiored  above,  the  details  of  the  mechanisms 


which  are  involved  are  quite  different,  and  revised  roughness  models  are  required 
in  order  to  include  these  types  of  surface  irregularities.  Again,  the  perturbed 
mean  flow  profiles  are  expected  to  lead  to  earlier  boundary-layer  transition. 
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Although  the  modification  of  the  boundary  layer's  mean  flow  profile  by  roughness 
elements  is  expected  to  be  the  dominant  effect,  the  increased  disturbance  level 
within  the  boundary  layer  is  also  important  and  should  be  included.  A parametric 
study  of  the  effect  on  the  boundary  layer  of  varying  intensity  and  spectral  con- 
tent of  the  disturbances  should  be  performed  and  the  connection  between  a given 
surface  condition  and  the  magnitude  of  these  parameters  needs  to  be  found. 

These  considerations,  however,  have  not  been  incorporated  into  the  present  in- 
vestigation for  lack  of  a definitive  physical  correspondence  between  the  surface 
conditions  and  the  factors  influencing  the  disturbance  level  within  the  boundary 
layer. 
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2.0  EFFECTS  OF  DISTRIBUTED  ROUGHNESS  ON  BOUNDARY-LAYER  TRANSITION 

2.1  A Brief  Description  of  the  Roughness  Model 

A physical  model  for  evaluating  the  effects  of  distributed  surface  roughness  on 
transition  characteristics  was  originally  developed  and  discussed  by  Kosecoff, 
Ko  and  Merkle  (1976).  In  order  to  facilitate  the  present  discussion,  a brief 
description  of  the  model  is  given  in  this  section.  Detailed  derivations  and 
the  background  of  this  particular  model  formulation  are  given  in  the  original 
reference. 


In  this  model,  the  roughness  elements  are  assumed  to  induce  an  enhanced  momentum 

and  heat  transfer  near  the  surface.  The  enhanced  momentum  and  heat  transfer  are 

incorporated  into  the  equations  of  motion  by  means  of  a momentum  diffusivity, 

e„,  and  a thermal  diffusivity,  The  momentum  diffusivity  is  assumed  to  be 
III  n 


^m  ” ^max 


e-ei(y/k)* 


(1) 


where  6i  is  a constant  (which,  effectively,  determines  the  "roughness  sub-layer 
thickness),  y is  the  distance  away  from  the  wall  in  the  normal  direction  measured 
from  the  bottom  of  the  roughness  elements,  k is  the  roughness  height,  and  is 
the  amplitude  of  the  momentum  diffusivity  at  the  wall  which  is  represented  by  the 
following: 


['  " “•’(-"V''*’]  ’ (2) 

Rej^  is  the  roughness  Reynolds  number  defined  by 

Rej^  = \k/V|^,  (3) 

where  the  subscript  k denotes  a quantity  evaluated  at  y = k,  K is  a constant, 

+ ^ + 
and  A is  a threshold  roughness  Reynolds  number.  Obviously,  for  , the 

effects  of  roughness  diminish  exponentially. 
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The  thermal  diffusivity  is  defined  in  terms  of  a turbulent  Prandtl  number,  Pr^, 
so  that 


t ■ 


(4) 
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The  two  quantities  and  are  incorporated  into  the  momentum  and  energy 
equations  by  means  of  an  effective  viscosity  and  an  effective  thermal  con- 
ductivity which  are  defined  as 


' M + pe„ 

(5) 

Kj  = K ♦ pCp 

(6) 

where  y and  K are  the  molecular  viscosity  and  thermal  conductivity,  p is  the  fluid 
density,  and  Cp  is  the  specific  heat.  To  incorporate  the  effects  of  surface  rough- 
ness in  the  boundary-layer  calculation,  the  viscosity,  y,  and  the  conductivity,  K, 
are  replaced  by  their  rough-wall  counterparts,  y^,  and  K^. 

In  the  numerical  evaluation  of  the  effect  of  distributed  roughness,  it  is  assumed 
that  the  turbulent  Prandtl  number  Is  approximately  one  (i.e.,  and  the 

following  values  have  been  determined  for  the  constants  appearing  in  the  model: 

= 0.094 

= 40  (7) 

Pi  = 1 

2.2  Validation  of  the  Model  - Achenbach’s  Experiment 

Very  little  experimental  information  on  the  effect  of  distributed  surface  rough- 
ness on  transition  has  been  reported.  It  is  generally  accepted  that  distributed 
roughness  causes  transition  to  occur  earlier  than  it  would  on  a smooth  wall,  but 
it  is  difficult  to  find  quantitative  measures  of  this  shift,  even  for  the  most 
simple  flow  fields.  Because  of  this  lack  of  data,  our  distributed  roughness 
model  has  been  compared  with  only  a limited  number  of  experimental  results, 
including  the  transition  data  reported  by  Feindt  (1957),  and  the  nosetip 
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transition  results  which  were  obtained  in  the  ABRES  reentry  program  (Merkle, 

1976).  In  the  comparisons  with  Feindt's  data,  the  model  provided  correct  quali- 
tative (as  well  as  reasonably  accurate  quantitative)  predictions  of  the  combined 
effects  of  roughness  and  pressure  gradient  (both  favorable  and  unfavorable)  on 
the  location  of  transition  (Feindt's  zero-pressure  gradient  results  were  used  to 
complete  the  calibration  of  the  constants  in  the  model).  With  regard  to  the 
reentry  vehicle  nosetip  data,  the  model  again  predicted  the  correct  qualitative, 
and  reasonably  accurate  quantitative,  effects  of  roughness  on  the  location  of 
transition  in  the  presence  of  several  simultaneous  effects  such  as  severe  surface 
cooling,  strong  pressure  gradients,  boundary- layer  blowing  at  the  surface  (simu- 
lated ablation),  and  changes  in  nosetip  geometry.  Additional  experimental  data 
which  have  not  previously  been  used  for  comparative  purposes  are  the  results 
reported  by  Achenbach  (1971).  Since  these  experimental  results  encompass  varia- 
tions of  nearly  two  orders  of  magnitude  in  the  unit  Reynolds  number,  they  repre- 
sent a meaningful  additional  test  of  the  ability  of  the  model  to  predict  the 
effects  of  distributed  roughness  on  the  location  of  transition.  Accordingly, 
the  model  has  been  applied  to  the  configuration  which  was  tested  by  Achenbach, 
and  the  predictions  of  the  model  have  been  compared  with  the  experimental 
measurements.  j 

The  Achenbach  Experiment 

Achenbach's  experiments  include  the  measurement  of  the  location  of  transition  on 
a circular  cylinder  immersed  in  a crossflow  at  different  freestream  Reynolds 
numbers.  The  Reynolds  number  based  on  the  diameter  of  the  cylinder  (150  mm)  was 
systematically  varied  between  6.0  x 10"  and  4.0  x 10®  so  that  the  location  of 
transition  ranged  from  less  than  10  degrees  from  the  forward  stagnation  line  to  I 

as  much  as  100  degrees  away  (on  the  leeward  side  of  the  cylinder).  Three  differ-  | 

ent  roughnesses  were  used  in  the  experiment;  the  first  two  were  constructed  by  H 

gluing  emery  paper  to  the  surface  of  the  cylinder,  while  the  third  was  constructed  H 

by  brazing  spherical  particles  to  the  surface.  The  roughness  heights  of  these  Q 

three  surfaces  were  estimated  by  generating  a turbulent  boundary  layer  over  the  I 

surface,  measuring  the  resulting  skin  friction,  and  equating  the  roughness  to  an  I 

equivalent  sand  grain  roughness  by  means  of  Nikuradse's  data.  The  roughness  I 
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heights,  as  detennined  in  this  manner,  were  k = 165,  675,  and  1350y  (ly=0.001nin) 
The  location  of  transition  was  detennined  by  measuring  the  local  surface  shear 
stress  by  means  of  a probe  imbedded  in  the  cylinder.  The  shear  stress  at  any 
angular  location  could  be  determined  by  rotating  the  cylinder  about  its  axis  so 
as  to  position  the  probe  at  the  desired  angular  location.  Since  the  flow  near  the 
stagnation  region  of  the  cylinder  represents  a very  strong  favorable  pressure 
gradient  flow  field  (starting  from  a Hartree  pressure  gradient  parameter  of  3=1.0), 
these  data  include  the  simultaneous  effects  of  roughness  and  pressure  gradient  on 
transition. 

Mean  Flow  Properties  of  the  Cylinder  Boundary  Layer 

All  transition  distances  quoted  in  Achenbach's  paper  are  normalized  by  the 
diameter  of  the  cylinder.  The  use  of  a local  boundary-layer  thickness  (such  as 
the  displacement  or  momentum  thickness)  should  give  considerably  more  insight  into 
the  experimental  data  because  these  boundary-layer  thicknesses  are  the  proper 
local  lengths  with  which  the  roughness  height  should  be  compared.  In  order  to 
estimate  the  boundary-layer  characteristics  for  Achenbach's  experiment,  the 
standard  formula  for  Incompressible,  inviscid  flow  about  a cylinder  was  used  to 
estimate  the  pressure  distribution.  The  boundary- layer  characteristics  were  then 
computed  from  a numerical  boundary-layer  code  which  was  a modification  of  the  code 
developed  by  Price  and  Harris  (1972).  Figure  1 through  4 show  some  of  the  charac- 
teristics of  the  boundary  layer  as  computed  for  the  smocth-wall  case.  A schematic 
of  the  flow  around  the  cylinder  is  shown  in  Figure  1 and,  as  shown  in  Figure  2, 
the  Reynolds  number  based  on  the  local  boundary-layer  displacement  thickness 
remains  quite  low  (generally  below  600)  in  the  region  where  transition  was 
observed.  It  is  interesting  to  note  that  these  Reynolds  numbers  are  smaller 
than  those  corresponding  to  transition  on  a smooth  flat  plate;  thus,  the 
presence  of  the  roughness  on  the  cylinder  more  than  offsets  the  effects  of  the 
strong  favorable  pressure  gradient. 

The  computed  momentum  thickness  of  the  smooth-wall  boundary  layer  is  shown  in 
Figure  3.  Note  that  the  momentum  thickness  is  generally  less  than  75vi,  while 
the  smallest  roughness  height  used  by  Achenbach  was  165y.  Thus,  the  roughness 
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tested  by  Achenbach  are  considered  quite  large.  The  third  of  the  smooth-wall 
calculations,  the  boundary-layer  shape  factor,  is  shown  In  Figure  4.  Note  that 
the  shape  factor  is  almost  constant  over  much  of  the  cylinder  for  the  smooth- 
wall  calculation. 


When  the  presence  of  the  165y  surface  roughness  Is  Included  In  the  boundary-layer 
computation,  the  shape  factor  Increases  considerably  and  varies  with  both  the 
freestream  Reynolds  number  and  the  angular  location  on  the  cylinder  as  Is  also 
shown  In  Figure  4.  The  roughness  Reynolds  number,  Rej^,  Is  shown  In  Figure  5. 
Again,  the  freestream  Reynolds  number  has  a strong  effect  on  the  roughness 
Reynolds  number.  Finally,  the  momentum  thickness  distribution  for  the  rough-wall 
boundary  layer  Is  shown  In  Figure  6.  By  comparison  with  Figure  3,  It  can  be  seen 
that  roughness  has  only  a small  effect  on  the  momentum  thickness  development,  so 
that  the  large  Increase  In  the  rough-wall  shape  factor  is  primarily  attributable 
to  the  change  In  the  displacement  thickness. 


Stability  Calculations  and  Transition  Predictions  for  Achenbach's  Experiments 


The  neutral  stability  curve  for  the  boundary-layer  on  Achenbach's  cylinder  with  a 
roughness  of  165y,  and  at  a Reynolds  number  of  Re^  = 4.0  x 10®  Is  presented  In 
Figure  7.  From  this  figure  It  may  be  shown  that  the  critical  Reynolds  number 
occurs  very  close  to  the  forward  stagnation  line  (at  about  ((>  * 5°).  When  the 
total  growth  rate  of  each  of  the  various  disturbance  frequencies  are  computed 
and  the  "e®"  transition  criterion  is  used,  the  transition  location  as  a function 
of  Reynolds  number  Is  obtained  for  the  165vi  case;  these  results  are  shown  In 
Figure  8 along  with  faired  curves  representing  the  experimental  data  for  the 
three  roughnesses  tested  by  Achenbach. 


The  comparison  shown  In  Figure  8 Indicates  that  the  present  roughness  model  over- 
predicts the  effects  of  roughness  on  transition  for  this  case.  Part  of  this 
overprediction  may  be  attributed  to  the  skin  friction  transition  criterion  which 
was  employed  by  Achenbach.  Skin  friction  Is  known  to  be  a relatively  late  Indi- 
cator of  transition  (that  Is,  the  Increase  In  skin  friction  generally  Is  observed 
toward  the  latter  stages  of  boundary-layer  transition),  whereas  the  "e*"  method 
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attempts  to  predict  the  starting,  or  early,  portion  of  the  transition  phenomenon. 

The  ambiguity  in  the  roughness  height  (owing  to  the  "roundabout"  procedure  which 
Achenbach  used  to  quantify  the  roughness  height)  also  introduces  an  uncertainty 
into  the  prediction.  For  a qualitative  assessment  of  the  trend  predicted  by  the 
present  model,  the  effects  of  these  uncertainties  could  be  lessened. 

Finally,  the  predicted  effect  of  changes  in  the  roughness  height  on  the  location 
of  transition  is  compared  with  the  experimental  measurements  in  Figure  9.  Again, 
the  analysis  predicts  the  correct  trend  for  the  effects  of  roughness  on  the 
transition  location,  but  it  quantitatively  overestimates  the  effects  of  the 
roughness. 

In  summary,  because  of  the  uncertainties  associated  with  the  Achenbach  experiment, 
a direct  quantitative  comparison  of  the  transition  location  could  not  be  made. 

The  calculations  do,  however,  demonstrate  a qualitative  agreement  with  the  experi- 
mentally observed  effect  of  surface  roughness  on  the  transition  of  a boundary 
layer  about  a circular  cylinder.  In  general,  the  predictions  given  by  the 
present  roughness  model  appear  to  be  somewhat  conservative. 

2.3  Application  of  the  Distributed  Roughness  Model  to  the  Boundary  Layer  on 
a Realistic  Body 

Some  computations  of  the  effect  of  distributed  surface  roughness  for  a realistic 
body  shape  have  been  completed  and  are  described  below.  A particular  vehicle  desig- 
nated as  Body  "H"  has  been  chosen  for  this  investigation.  The  body  contour  results 
in  a fairly  strong  favorable  pressure  gradient  over  most  of  its  length;  details  of 
the  body  geometry  and  the  (inviscid)  surface  velocity  are  given  in  Figure  10. 

For  this  body  it  was  assumed  that  the  surface  temperature  was  11°C  above  the 
ambient  water  (i.e.,  AT  = 11®C)  and  that  the  Reynolds  number  based  on  the  overall 
length  was  300  million.  Different  roughness  levels  were  investigated  to  find  a 
‘'threshold"  level  at  which  the  roughness  began  to  have  a non-negl igible  effect  on 
the  stability  characteristics  of  the  boundary  layer.  The  streamwise  variation  in 
the  Reynolds  number  based  on  displacefnent  thickness  is  shown  in  Figure  11  for  both 
the  smooth-wall  case  and  for  a surface  roughness  of  25y.  As  can  be  seen,  surface 
roughness  causes  a very  modest  increase  in  6*  compared  to  the  smooth-wall  case. 
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The  momentum  thickness  development  for  this  body  is  shown  in  Figure  12,  and  the 
ratio  of  the  roughness  height  to  momentum  thickness  is  shown  in  Figure  13.  As 
can  be  seen  from  this  latter  figure,  the  growth  of  the  boundary  layer  from 
its  relatively  thin  value  at  the  stagnation  point  indicates  that  a threshold 
roughness  level  's  a function  of  distance  along  the  body.  Small  roughness  can 
have  substantial  effects  near  the  forward  portion  of  the  body  but  they  will  be 
effectively  buried  beneath  the  boundary  layer  by  the  middle  of  the  body.  In 
particular,  note  that  for  the  25p  surface  roughness  case  k/0>  0.2  (the  value 
quoted  as  a roughness  threshold  in  the  similar-flow  investigations  of  our  pre- 
vious report)  only  over  the  first  15  percent  of  the  body  length;  the  12. 5p  sur- 
face roughness  results  in  k/e>  0.2  only  over  the  first  2.5  percent  of  the  body 
length;  and  the  2.5p  surface  roughness  case  never  results  in  k/0>  0.2.  The 
corresponding  roughness  Reynolds  numbers  are  shown  in  Figure  14.  As  can  be  seen 
from  this  figure,  the  roughness  Reynolds  number  increases  nearly  with  the  square 
of  the  roughness  height  so  that  small  changes  in  the  roughness  height  can  result 
in  substantial  changes  in  the  roughness  Reynolds  number.  Accordingly,  small 
changes  in  the  roughness  height  can  have  a sizeable  effect  on  the  location  of 
transition.  Roughness  Reynolds  numbers  above  30  are  noted  for  the  25y  surface 
roughness  case. 

The  stability  characteristics  for  the  boundary  layer  about  this  body  for  several 
roughness  heights  are  shown  in  Figures  15  through  18.  Figure  15  shows  the  neutral 
stability  curve  for  the  smooth-wall  case.  When  a roughness  height  of  16p  is 
assumed,  the  unstable  region  distorts  slightly  as  shown  in  Figure  16.  A compari- 
son between  Figures  15  and  16  indicates  that  roughness  has  caused  the  critical 
Reynolds  number  to  occur  closer  to  the  stagnation  point,  and  has  increased  the 
extent  of  the  unstable  region,  but  that  the  downstream  portions  of  the  unstable 
region  are  virtually  unchanged.  This  insensitivity  to  the  roughness  at  the 
farther  downstream  stations  occurs  because  the  roughness  becomes  embedded  in  the 
boundary  layer  as  it  thickens. 

When  a surface  roughness  of  25y  is  assumed,  the  critical  Reynolds  number  occurs 
nearly  at  the  stagnation  point  and  a large,  highly-unstable  region  exists  as  shown 
in  Figure  17.  A comparison  between  Figures  16  and  17  indicates  the  strong  effect 
of  increasing  roughness  height  on  the  boundary-layer  stability  characteristics. 
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Finally,  the  integrated  disturbance  growth  rates  for  these  two  cases,  the  16u  and 
the  25u  surface  roughnesses,  are  shown  in  Figure  18.  As  can  be  seen,  the  16u 
surface  roughness  gives  rise  to  disturbance  amplifications  as  large  as  e®.  The 
25u  surface  roughness,  however,  yields  an  amplification  in  excess  of  e“  for  all 
three  disturbance  frequencies  shown.  This  last  figure  clearly  indicates  that  the 
threshold  roughness  for  this  body  operating  at  this  speed  (i.e.,  the  Reynolds 
number  u L/v  = 300  x 10®)  and  with  a wall  overheat  of  11°C  is  about  IGp. 
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3.0  INCORPORATION  OF  THE  ROUGHNESS  MODEL  INTO  THE  TAPS  CODE 

As  described  in  the  previous  section,  the  model  for  estimating  the  effects  of 
distributed  surface  roughness  on  boundary-layer  transition  has  been  developed  and 
validated  with  experimental  data  with  qualitatively  favorable  results.  The  model, 
however,  needs  to  be  validated  further  with  a few  well -control led  quantitative 
measurements  before  confidence  in  the  model  and  its  associated  "constants"  can  be 
fully  established.  At  the  present  time,  the  model  may  be  used  as  a qualitative 
(possibly  conservative)  guide  for  estimating  the  effects  of  distributed  surface 
roughness  for  design  purpose.  Therefore,  in  order  to  make  the  roughness  model 
more  readily  available  to  the  other  investigators  in  the  laminar  flow  research 
area,  an  effort  has  been  devoted  in  this  study  to  incorporate  the  roughness  - 
model  into  the  TAPS  code. 


3.1  Analytical  Approach 

The  boundary- layer  segment  of  the  TAPS  program  is  based  on  the  Cebeci -Smith  (1974) 
finite-difference  boundary-layer  program.  The  general  solution  procedure  used  in 
the  program  is  shown  in  the  diagram  below. 


Repeat  this  cycle  until 

then  proceed  to  the  next  station. 
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where  f"  is  the  velocity  gradient  at  the  wall,  i is  the  iteration  number  and  a is 
w 

a convergence  criterion. 

As  mentioned  in  Subsection  2.1,  the  effects  of  distributed  surface  roughness  are 
incorporated  by  replacing  the  molecular  viscosity  and  thermal  conductivity,  p and 
K,  by  their  rough-wall  counterparts,  and  in  the  momentum  and  energy  equations. 
Hence,  the  normal iaed  boundary-layer  equations  for  either  two-dimensional  or 
axisymmetric  flow  can  be  written  as  follows  (Gentry,  1976): 


Continuity 


^ (r"p.)  ♦ ^ (rV)  ■ 0 


Momentum 


'I  ® 


Energy 


where 


4.  3H  . 

P“  97  * P’  17  • 


u streamwise  (x)  component  of  fluid  velocity 
V cross-stream  (y)  component  of  fluid  velocity 
r radial  (or  normal)  distance  from  the  axis 
p density 

n = 0 for  two-dimensional  flow  and 
= 1 for  axisymmetric  flow 

H total  enthalpy 

p^  total  or  effective  viscosity 

K.  total  or  effective  thermal  conductivity 

Cp  specific  heat  at  constant  pressure 

( )g  conditions  at  the  outer  edge  of  the  boundary  layer 


t i ! 


In  order  to  remove  the  singularity  at  x*0,  and  to  stretch  the  coordinate  normal  to 
the  flow  direction,  the  Levy-Lees  (Hayes  and  Probstein,  1959)  and  the  Probsteln-Elllot 
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transformations  (1956)  are  used  in  obtaining  the  solution  of  the  equations.  The 
transformations  are 

“ Pe^e^e  T 


pu  i n 

dn  = ^ dy 

(25)’*  ' '■ 


r = rQ  + y cos  a (13; 

where  r^  is  the  body  radius  and  a is  the  angle  that  the  surface  makes  with  the 
body  axis. 

To  simplify  the  procedure  for  incorporating  the  roughness  model,  the  roughness 
height  k is  expressed  in  the  transformed  coordinate  n.  Thus,  at  a given  stream- 
wise  position,  the  value  of  n at  the  top  of  the  roughness  is 


' 1(25)**  i 


Pe“e  r k^ 

c U + f cos 

L(25)^  L 0 2 


J L(2C)^  I-  ‘’-I 


Once  the  roughness  height  is  defined  in  the  transformed  coordinate  system,  its 
effect  can  be  computed  at  each  streamwise  location.  The  procedure  used  for 
incorporating  the  roughness  model  into  the  TAPS  program  is  described  below  in 
conjunction  with  Figure  19. 

1.  The  computation  will  include  the  effect  of  surface  roughness  If  the 
control  flag  LG23  = 3 in  the  input  data. 

2.  Input  data  should  include  the  roughness  height  k along  the  body  surface. 

3.  Calculate  the  roughness  height  k in  the  transformed  coordinate  n> 
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4.  Obtain  the  variations  of  velocity  u,  temperature  T,  viscosity  y, 
density  p in  the  cross-stream  (n)  direction  from  the  TAPS  boundary- 
layer  flow  calculation. 

5.  Calculate  the  roughness  Reynolds  number  Rej^. 

6.  Calculate  the  momentum  diffusivity  and  thermal  diffusivity 
according  to  the  roughness  model. 

7.  Obtain  the  total  or  effective  viscosity  y^  and  thermal  conductivity  Kj. 

8.  Replace  the  molecular  viscosity,  y,  and  conductivity,  K,  by  their  rough- 
wall  counterparts  y^  and  ICj.  in  the  momentum  and  energy  equations. 

9.  Resume  the  boundary-layer  calculation  in  the  TAPS  computer  program. 

3.2  Descriptions  of  the  Changes  in  the  TAPS  Code 

Modification  of  the  Input  Data 

The  Input  Data  format  is  the  same  as  the  original  TAPS  Input  format  except  for 
the  following  changes: 

i)  Flag  Control  Card 

Column  Code  Routine  Format  Explanation 

23  LG23  INP2  II  Laminar  flow  roughness  control  flag 

= 0 smooth  surface 

= 3 roughness  height  k is  input  on  the 
boundary-layer  station  data  card 

ii)  Boundary-Layer  Station  Data  Cards 

Column  Code  Routine  Format  Explanation 

41-46  RP(I)  INP2  F6.0  Equivalent  sand-grain  roughness  height 

for  laminar  flow,  inches 

(NOTE:  The  original  TAPS  Code  has  the  same  input  format  for  turbulent  flow 
calculations.) 

A sample  of  the  input  data  is  shown  in  Table  1. 


Modification  of  the  TAPS  Program 


The  structure  of  the  TAPS  Boundary-Layer  Program  is  presented  in  Table  2.  The 
primary  modifications  to  the  computer  program  are  in  the  FLP2  Subroutine  for  the 
calculation  of  the  fluid  properties.  These  modifications  are  shown  in  Tables  3-1 
and  3-2.  Subroutine  OTPT  also  has  been  modified  as  shown  in  Table  3-3  in  order 
to  simplify  the  output  format.  A sample  of  the  output  summary  is  presented  in 
Table  4. 


3.3  Discussion  of  the  Results 


Calculations  on  the  effects  of  distributed  roughness  for  the  Body  "H"  have  been 
performed  by  using  the  modified  TAPS  code.  A comparison  of  the  results  with  those 
calculation  in  Section  2 is  presented  in  Figure  20  through  Figure  22. 

In  general,  results  from  the  two  codes  are  in  fairly  good  agreement.  For  the 
smooth-wall  case,  the  variation  in  the  two  computed  momentum  thicknesses  is  only 
2%  at  x/L  * 0.4.  The  effects  of  surface  roughness  on  the  roughness  Reynolds 
p'miber  Rej^  and  the  parameter  k/0  are  presented  in  Figures  21  and  22.  For  the 
small  roughness  height  k * 12. 5y,  the  discrepancies  from  the  two  codes  are  neg- 
ligible. For  the  roughness  height  k = 25v,  the  resulting  Rej^  and  k/0  from  Harris' 
code  are  slightly  higher  than  the  results  from  TAPS.  In  Subsection  2.2,  it  was 
concluded  that  the  analysis  from  Harris'  code  quantitatively  overpredicted  the 
effects  of  roughness  at  k = 165y  (the  experimental  data  of  Achenbach).  There- 
fore, by  intuition  it  is  anticipated  that  the  roughness  model  in  the  TAPS  code 
should  yield  a better  agreement  with  experimental  data  than  Harris'  code. 
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4.0  SOME  PRELIMINARY  CONSIDERATIONS  ON  OTHER  TYPES  OF  SURFACE  ROUGHNESS 


Thus  far,  the  cnalytical  effort  at  assessing  the  effects  of  surface  roughness  on 
the  location  (ff  boundary- layer  transition  has  dealt  with  distributed  roughness 
only.  The  basic  mechanism  envisioned  in  that  model  is  an  enhanced  momentum  and 
heat  transfer  near  the  wall  as  a result  of  the  unsteady  flow  downstream  of  each 
distributed  roughness  element.  This  effect  of  surface  roughness  is  manifested 
through  the  change  of  the  stability  characteristics  of  the  modified  mean  flow 
profiles  within  the  laminar  portion  of  the  boundary  layer.  An  implicit  require- 
ment of  this  model  is  the  presence  of  a separated  flow  region  behind  each  rough- 
ness element.  A lower  bound  cut-off  of  Rej^  = >J|^k/V|^  = 40  has  been  used  in  the 
model  to  account  for  the  non-existence  of  a separated  region,  and,  therefore, 
the  absence  of  unsteady  flow  at  roughness  Reynolds  numbers  below  this  value. 

This  characteristic  value  was  assumed  to  be  the  same  as  the  incipient  Reynolds 
number  for  the  unsteady  wake  flow  behind  a cylinder.  The  arbitrariness  of  this 
cut-off  and  the  need  for  examining  the  effects  of  surface  roughness  without  flow 
separation  in  order  to  assess  the  minimum  critical  roughness  height  with  more 
confidence  have  led  us  to  examine  and  to  formulate  the  problem  for  the  effect  of 
a wavy  wall  as  summarized  below.  The  intent  is  to  provide  a lower  bound  on  the 
effect  of  roughness  as  well  as  to  simulate  the  frequently  postulated  machining- 
induced,  organized  surface  roughness. 

For  larger  roughness  heights  under  practical  conditions,  a model  of  isolated  two- 
dimensional  and/or  three-dimensional  roughness  elements  is  more  representative. 
Their  effect  on  boundary-layer  transition  is  often  observed  to  be  significant 
and  to  dominate.  In  Subsection  4.2,  we  present  a brief  review  and  attempt  to 
investigate  the  effect  of  an  isolated  two-dimensional  roughness  element  using 
the  triple-deck  formulation  of  Smith  (1973). 


Recently,  Lessen  and  Gangwani  (1976)  have  presented  an  analysis  of  the  flow  over 
a wavy  wall.  Their  analysis  restricts  the  wavelength  of  the  roughness  to  be  small 
compared  with  the  wavelength  of  a Tollmien-Schlichting  wave,  and,  in  addition, 
restricts  the  height  of  the  roughness  so  that  the  wall  angles  are  shallow  and 
separation  does  not  occur  behind  the  individual  bump.  Viscous  flow  over  a wavy 
wall  causes  a mean  distortion  to  the  flow  field  (if  second  order  tenns  are  kept) 
which  will,  in  turn,  affect  the  stability  characteristics  of  the  boundary  layer. 
Since  this  mechanism  for  distorting  the  mean  flow  is  a steady-flow  mechanism,  its 
effect  is  not  included  in  the  distributed  roughness  model  which  applied  to  rough- 
ness elements  whose  height  is  comparable  to  their  characteristic  wavelength.  The 
unsteady  flow  effects  of  the  distributed  roughness  model  are  expected  to  dominate 
the  second-order  steady  flow  distortions  which  are  predicted  by  Lessen' s analysis 
when  the  roughness  is  large,  but,  in  the  limit  when  the  roughness  heights  are 
small,  the  latter  effect  may  be  the  more  important.  The  purpose  of  reviewing 
Lessen's  analysis  is  to  determine  whether  his  analysis  predicts  a more  signifi- 
cant effect  of  roughness  in  the  very  small  roughness-height  region  where  laminar 
flow  underwater  bodies  are  expected  to  operate. 

In  the  wavy-wall  analysis,  the  flow  field  is  divided  into  three  separate  parts, 
a distorted,  spatially  averaged  mean  flow,  a steady  perturbation  at  a wavelength 
corresponding  to  that  of  the  wavy  wall,  and  an  unsteady  Tollmien-Schlichting  wave. 
In  order  to  keep  the  analysis  simple,  and  to  cast  emphasis  upon  the  effects  of 
the  wall  waviness  on  the  stability  characteristics  of  the  flow,  the  T-S  waves 
will  be  treated  as  being  infinitesimal  in  strength  (i.e.,  they  are  considered 
as  being  linear  in  a traditional  "e”'  sense),  while  the  roughness  amplitude  is 
taken  as  being  small,  but  finite.  In  the  analysis,  the  distorted  mean  flow  will 
be  denoted  by  a variable  with  an  over-bar,  the  periodic  perturbation  due  to  the 
roughness  will  be  denoted  by  a superscript  asterisk,  and  the  T-S  wave  will  be 
indicated  by  a prime.  Other  symbols  are  taken  as  standard:  u and  v are  the 
velocity  components  in  an  x-y  cartesian  coordinate  system,  the  pressure  divided 


by  the  density  1s  Indicated  by  p,  and  the  kinematic  viscosity  by  v.  Thus,  the 
streamwlse  velocity  component  Is  written  as 


u = u + u*  + u' 


(15) 


If  expressions  of  the  form  of  Equation  (15)  are  Introduced  Into  the  Navler-Stokes 
equations,  and  the  flow  field  Is  averaged  over  time  and  over  one  wavelength  of 
the  roughness,  the  equations  for  the  mean  flow  become 
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where  the  boundary-layer  approximations  have  been  applied,  and  the  second  order 
quantities  In  the  roughness-induced  distortion  (but  not  the  T-S  wave  as  Indicated 
above)  have  been  retained.  The  symbol, <>,  Is  used  to  Indicate  the  spatial  average 
over  one  wavelength  of  the  roughness.  By  again  returning  to  the  complete  Navler- 
Stokes  equations.  Introducing  Equation  (15),  time  averaging,  subtracting  Equations 
(16)  and  (17),  and  applying  parallel-flow  approximations,  we  obtain  the  equations 
governing  the  perturbation  Induced  by  the  roughness: 


= If  * f * I?  ■ 
= lf 


(18) 

(19) 
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(20) 


This  system  of  equation  Is  Identical  to  that  which  governs  Tollmlen-Schllchting 
waves,  except  that  the  unsteady  terms  are  not  present.  As  will  be  seen  later, 
the  boundary  conditions  corresponding  to  this  equation  system  are  non-homogeneous 
so  that  non-trivlal  solutions  to  Equations  (18)  and  (20)  can  be  found.  Note  that 
a linear  approximation  tq  this  set  of  equations  Is  sufficient  In  order  to  main- 
tain the  Reynolds  stress  term  which  appears  In  the  mean  flow  equation  (16);  the 
Reynolds  stress  term  <u*v*>.  Is  non-zero  because  of  the  presence  of  the  mean 
shear  u (y). 
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The  equations  for  the  Tonmien-Schllchting  wave  are  obtained  by  standard 
techniques  and  have  the  standard  form  except  that  the  mean  flow  coefficients 
which  appear  in  the  Orr-Sommerfel d equation  are  those  of  the  distorted  mean 
flow  as  computed  from  Equations  (16)  and  (17). 

The  analysis  of  the  flow  over  a wavy  wall  is  completed  by  specifying  the  boundary 
conditions  for  the  three  systems  of  equations.  These  boundary  conditions,  which 
can  be  obtained  by  expanding  the  solution  in  a Taylor's  series  at  the  wall,  are 
as  follows:  the  boundary  condition  for  the  mean  flow  is 

u(0)  = - j <y^>  Uyy(O)  - <y^^uj(0)>  (21) 

where  y^  is  the  amplitude  of  the  roughness;  the  boundary  conditions  for  the  rough- 
ness induced  disturbance  is 

u*(0)  = y„Uy(0)  ; (22) 

while  the  boundary  conditions  for  the  Tollmien-Schlichting  wave  are  the  standard 
homogeneous  ones,  both  at  infinity  and  at  the  wavy  wall  surface. 

Some  effects  of  the  wavy  wall  on  the  critical  Reynolds  number  are  given  in 
Figure  23.  These  results  have  been  computed  by  Lessen;  results  for  higher  Rey- 
nolds numbers  are  not  available.  As  can  be  seen,  roughness  is  destabilizing, 
and,  if  the  limits  of  the  analysis  are  severely  strained,  changes  of  about  10 
percent  can  be  noted  in  the  critical  Reynolds  number. 

Lessen 's  model  has  been  extended  to  include  tbe  effects  of  pressure  gradient  and 
surface  heat  transfer.  Some  preliminary  order  of  magnitude  estimates  have  been 
made  for  the  effect  of  surface  heating  combined  with  a wavy  wall.  These  results, 
which  are  for  wavelength-to-height  ratios  which  ensure  that  the  flow  remains 
attached,  have  indicated  that  heating  has  a very  weak  effect  on  stability  and 
transition  in  this  circumstance.  Consequently,  our  net  assessment  for  the  effect 
of  surface  roughness  on  stability  is  that  even  in  the  limit  of  very  small 
"roughness"  heights,  waviness  never  dominates  the  distributed  roughness  effect 
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(for  similar  values  of  roughness  height).  Therefore,  we  may  consider  the  threshold 
roughness  height  established  with  the  distributed  roughness  model  to  be  conserva- 
tive. Further  investigation  of  the  effect  of  a wavy  wall  appears  unwarranted  at 
the  present  time. 

4.2  Single  Two-Dimensional  Bump 

The  effects  of  a single  two-dimensional  bump  on  the  boundary-layer  mean-flow 
profile  have  been  considered  by  Smith  (1973)  who  used  a "triple-deck"  analysis 
to  determine  its  effects.  The  triple  deck-analysis  expands  the  equations  in 
terms  of  a small  parameter  e,  which  is  defined  as 

e = Re'^  (23) 

where  6 is  the  local  boundary-layer  thickness.  The  height  of  the  bump  is  of  order 
£6  and  its  width  must  be  of  order  6/e.  Asymptotic  equations  are  written  for  each 
of  these  regions  and  proper  matching  procedures  are  used  to  determine  the  inter- 
action between  the  three  regions. 

In  Smith's  analysis,  the  inner  region  extends  from  the  wall  to  a height  which  is 
similar  in  order  of  magnitude  to  that  of  the  bump,  0(e6).  The  main  or  middle 
deck  constitutes  the  remainder  of  the  boundary  layer  and  has  a thickness  of  the 
order  of  6.  The  upper  deck  extends  to  a height  of  the  order  of  6/e,  and,  being 
outside  the  boundary  layer,  accounts  for  the  interaction  between  the  inviscid 
flow  and  the  bump. 

The  detailed  results  of  the  triple-deck  analysis  show  that  the  predominant 
effect  of  the  bump  on  the  boundary  layer  is  to  cause  a local  distortion  in  the 
impressed  pressure  gradient.  The  effects  of  the  bump  are,  however,  largely  local 
and  decay  to  small  levels  within  a length  on  the  order  of  a bump  width.  The  net 
distortion  in  the  profile  which  exists  downstream  of  the  bump  is,  therefore, 
quite  small.  The  analysis  does,  however,  show  that  the  boundary-layer  formula- 
tion is  adequate  for  flows  over  bumps  of  this  nature,  but  that  the  interaction 
with  the  inviscid  flow  must  be  included  in  the  analysis.  The  details  of  includ- 
ing the  viscous/inviscid  interaction  require  a level  of  difficulty  which  is 


-21- 


almost  identical  to  computing  the  complete  boundary-layer  equations  numerically. 
Consequently,  the  effects  of  the  bump  can  equally  well  be  included  by  using  a 
standard  numerical  boundary- layer  scheme  and  using  proper  step  size  resolution 
in  the  vicinity  of  the  bump.  In  so  doing,  the  perturbation  to  the  local  pres- 
sure gradient  must,  of  course,  be  included  in  the  computation. 

In  the  recent  Westinghouse  tow-tank  tests  the  effect  of  a single  two-dimensional 
bump  such  as  that  described  above  was  examined.  The  following  is  an  analysis  of 
its  effect  on  the  boi.ndary  layer  and  its  transitional  characteristics. 

Effects  of  a Single  Two-Dimensional  Bump  on  Transition 

The  effects  of  a single  bump  on  the  location  of  transition  may  be  handled  quite 
differently  from  an  analytical  viewpoint  than  the  case  of  multiple  bumps  (dis- 
tributed roughness  or  wall  waviness).  For  the  single  bump  case,  it  is  possible 
to  resolve  the  details  of  the  distortions  near  the  location  of  the  bump  and 
include  them  in  both  the  boundary  layer  and  stability  calculations.  Consequently, 
as  an  augmentation  of  our  wavy-wall  and  distributed  roughness  analyses,  we  are 
also  reviewing  the  effects  of  a single  two-dimensional  bump  on  the  location  of 
transition. 

For  this  analysis,  rather  than  consider  a general  case,  we  have  restricted  our 
attention  to  a specific  bump  of  given  geometry  on  a particular  body.  The 
specific  bump  which  we  are  considering  is  the  one  which  was  installed  upon  the 
Reichardt  body  in  the  recent  MBT*  tow  tank  tests  (Westinghouse,  1977).  Our 
intent  is  to  compare  the  analytical  predictions  of  its  effect  on  transition 
with  the  observed  experimental  results.  The  geometry  of  this  particular  bump 
and  its  location  on  the  vehicle  are  shown  in  Figure  24.  As  can  be  seen  from  the 
figure,  the  height  of  the  bump  is  0.125  mm  and  its  length  is  about  38  mm.  At 
the  critical  velocity  of  5.69  m/sec  the  undisturbed  boundary- layer  momentum 
thickness  at  this  location  is  computed  to  be  about  0.15  mm.  Thus,  the  bump 
height  is  slightly  smaller  than  the  local  momentum  thickness  of  the  boundary 
layer,  and  is  about  250  boundary-layer  (momentum)  thicknesses  long. 


A single  bump  of  this  character  can  affect  the  location  of  transition  through  at 
least  three  distinct  mechanisms.  First  of  all,  the  bump  introduces  a local  per- 
turbation to  the  boundary- layer  velocity  field  which  alters  the  stability  charac- 
teristics of  the  boundary  layer.  This  aspect  of  a bump  was  reviewed  In  terms  of 
a triple-deck  analysis  (Smith,  1P73)  which  Indicated  that  the  use  of  boundary- 
layer  theory  In  the  vicinity  of  the  bump  was  acceptable  so  long  as  the  local 
pressure  perturbation  Is  taken  into  account.  The  second  effect  Is  that  the 
bump,  which  was  composed  of  silver  epoxy,  may  not  have  reached  the  same  surface 
temperature  overheat  as  the  rest  of  the  wall  even  though  It  Is  a reasonably  good 
heat  conductor.  Finally,  the  local  geometrical  characteristics  of  the  surface 
which  are  Introduced  by  the  bump  have  regions  of  concave  curvature  near  the 
front  and  the  back  sides  of  the  bump,  and  these  regions  could  stimulate  the 
growth  of  Taylor-Gortler  vortices  which.  In  turn,  could  lead  to  transition  at 
substantially  lower  amplification  ratios  than  the  "e®"  criterion. 


The  profile  of  the  surface  wave  which  Is  shown  in  Figure  24  can,  to  a first 
approximation,  be  represented  as 


y = I (1  - cos  ^S) 


(24) 


Z 


where  the  bump  height,  h.  Is  0.125  mm  and  the  bump  length,  X,  Is  38.0  mm.  Con- 
sidering first  the  Taylor-Gortler  instability,  we  compute  a maximum  concave 
curvature  of  c = 1.72  x 10®  mm~*  at  S = 0 and  X.  For  a vehicle  velocity  of 
5.69  m/sec  (11  knots),  the  Gortler  parameter,  Gj.  (=  Re0V^),  at  this  location 
Is  G^  * 12.7.  Stability  calculations  Indicate  Taylor-Gortler  Instabilities 
become  manifest  when  the  parameter  G^  exceeds  a value  of  1 or  2,  and  this  would 
suggest  that  Gortler  Instabilities  could  be  present  near  the  bump.  We  note, 
however,  that  the  sinusoidal  representation  of  the  shape  of  the  bump  Is  only  an 
approximation  and  that  the  actual  shape  (and,  hence,  the  curvature)  remains 
undetermined.  Nevertheless,  the  above  calculation  should  have  the  correct  order 
of  magnitude,  and  a strong  likelihood  of  Instability  Is  Indicated.  The  reason 
for  the  possible  appearance  of  a Taylor-Gortler  Instability  Is,  of  course, 
directly  related  to  the  high  value  of  Reg  for  the  boundary  layer  which  was 
computed  for  the  body  flow;  this  sensitivity  Is  conmon  In  all  high  Reynolds 
number  tests. 
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’ I ^ The  pressure  coefficient  in  the  neighborhood  of  the  bump  is  shown  in  Figuie  25. 

In  the  absence  of  the  bump,  the  pressure  coefficient  is  virtually  constant,  but 
with  the  bump  present,  a characteristic  adverse  pressure  gradient,  followed  by 
a favorable  and  then  a second  adverse  pressure  gradient  region  results.  When 

j translated  into  boundary-layer  characteristics,  the  boundary  layer  is  actually 

<5 

I thinned  somewhat  in  the  favorable  pressure  gradient  region,  as  shown  on  Figure  26, 

but  the  net  effect  of  the  bump  is  to  leave  the  boundary  layer  somewhat  thicker 

< do’vnstream  of  the  bump  as  compared  to  the  undistorted  case. 

O 

Some  preliminary  stability  calculations  in  which  the  pressure  perturbation  has 
been  included,  and  for  which  it  is  assumed  that  no  surface  heat  transfer  takes 
place  in  the  vicinity  of  the  bump,  have  been  obtained.  These  results  indicate 
^ that  the  effect  of  the  two  perturbations  taken  together  (no  surface  heating  plus 

pressure  gradient  effects)  are  sufficient  to  cause  a noticeable  change  in  the 
stability  characteristics  and  to  force  transition  to  occur  at  the  bump.  These 
preliminary  analyses  should  be  further  expanded  in  the  future  in  order  to  obtain 
a more  definitive  assessment  of  the  effects  of  an  isolated  bump  on  boundary- 
layer  transition. 


O 


5.0  SUMMARY 


In  this  report,  we  have  presented  a summary  of  our  recent  investigations  of  the 
effects  of  surface  roughness  on  boundary -layer  transition.  The  primary  effort 
has  focused  on  the  formulation,  validation  and  application  of  the  enhanced 
momentum  and  heat  transfer  model  for  distributed  roughness  elements.  Without 
further  quantitative  experimental  confirmation  of  this  model,  it  is  our  present 
opinion  that  this  model,  which  has  been  incorporated  into  the  TAPS  code,  can  pro- 
vide a reasonable  prediction  of  the  effects  of  randomly  distributed  surface 
roughness  on  the  transition  location  for  a realistic  underwater  vehicle.  The 
model,  in  particular,  points  out  the  reduced  effectiveness  of  surface  heating 
for  delaying  boundary-layer  transition  in  the  presence  of  surface  roughness. 

The  second  type  of  surface  irregularity  considered  was  waviness  which  could 
result  during  vehicle  construction.  The  model  of  Lessen  and  Gangwani  has  been 
extended  to  provide  a qualitative  assessment  of  this  effect.  The  preliminary 
results,  however,  suggest  a weak  influence  for  expected  surface  waviness  param- 
eters such  as  when  the  wave  height  to  wave  length  ratio  is  small.  Additional 
effort  is  needed  to  examine  the  effect  of  surface  waviness  in  a more  quantitative 
manner  and  should  include  relatively  high  wave  height  to  wave  length  ratios  for 
which  the  flow  may  no  longer  remain  attached  everywhere  on  the  surface. 

The  third  type  of  surface  irregularity  considered  was  a single,  two-dimensional 
bump.  The  triple-deck  formulation  of  Smith  was  used  to  investigate  the  effect 
on  the  mean  flow  profiles  within  the  boundary  layer  and,  therefore,  the 
stability  characteristics.  This  analysis  is  quite  preliminary  in  nature  and 
additional  effort  is  needed  before  any  quantitative  assessment  can  be  made. 

In  suitmary,  we  believe  that  surface  irregularities  can  play  a dominant  role  in 
determining  the  practical  range  of  applicability  of  advanced  laminar- flow  tech- 
niques. Extreme  care  should  be  taken  in  the  interpretation  of  experimental 
results  if  quantitative  information  on  the  surface  finish  conditions  is  not 
known.  The  present  phenomenological  model  for  distributed  surface  roughness, 
however,  can  be  used  for  a first-order  estimate  of  the  transition  location. 


Based  on  the  qualitative  comparison  with  the  existing  experimental  results,  we 
would  deduce  that  the  present  model  would  provide  a rather  conservative  estimate 
of  the  effects  of  surface  roughness.  A more  quantitative  assessment,  however, 
must  rely  on  additional  validations  with  carefully  performed  surface  roughness 
experiments.  These  experiments  are  critically  needed  for  an  improved  under- 
standing of  the  effects  of  surface  roughness  on  the  boundary-layer  transition. 
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’ f Table  2.  Structure  of  the  TAPS  Boundary-Layer  Program 
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Figure  5.  Roughness  Reynolds  Number  on  Cylinder 
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Figure  6.  Momentum  Thickness  on  Rouah  Cylinder 
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Figure  10.  Geometry  and  Velocity  Distribution  for  Body 


Effect  of  Roughness  on  Displacement  Thickness 
Reynolds  Number 


Variation  of  Roughness  Height  to  Momentum  Thickness  Ratio 


Variation 
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Disturbance  Amplification  Ratio  Within  Boundary  Layer 
on  Heated  Body  with  Rough  Wall 


eJC.n)  - 0.094  ^ Re^[l  - Exp  (-Re^/40)]  Exp  [-  (n/n,^)*] 
£^(5.11)  ■ 0.094  ^Re|^[l  - Exp  (-Re|^/40)]  Exp[-  (n/n,^)^] 


Figure  19.  Flow  Diagram  for  Incorporating  Roughness 
Model  into  the  TAPS  Code 


Figure  20.  Comparison  of  Momentum  Thickness 


of  Roughness  Reynolds  Number 
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Figure  23.  Effect  of  Wavy  Wall  on  Critical  Reynolds  Number 
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Figure  26.  Effect  of  Single  Bump  on  the  Boundary- Layer 

Displacement  Thickness  and  the  Momentum  Thickness 
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